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Abstract 

The time evolution of a two-level atom which is simultaneously exposed to the field of a 
running laser wave and a homogeneous gravitational field is studied. The result of the coupled 
dynamics of internal transitions and center-of-mass motion is worked out exactly. Neglecting 
spontaneous emission and performing the rotating wave approximation we derive the complete 
time evolution operator in an algebraical way by using commutation relations. The result is 
discussed with respect to the physical implications. In particular the long time and short 
time behaviour is physically analyzed in detail. The breakdown of the Magnus perturbation 
expansion is shown. 



PACS: 42.50.Vk 



1 Introduction 

Over the last ten years the manipulation of neutral atoms by laser light was dramatically improved and has 
led to the new field of atom optics (see Ref. [jjj and references therein). It includes in particular atomic 
intcrferometry and laser cooling of atoms, i.e., the preparation of a dense cloud of atoms with a narrow 
momentum distribution. The width of the atomic velocity distribution can be made as small as 1 cm/s 
This is of interest for atomic intcrferometry because the possibility to use slow atom beams enlarges the 
phase shifts caused by a broad class of external potentials ||. It is obvious that for atoms moving with 
a velocity of a few centimeters or meters per second for a time period of several milliseconds or more the 
influence of the earth's acceleration becomes important and cannot be neglected. For this reason it is of 
interest to study the time evolution of an atom moving in the gravitational field of the earth and the field 
of a laser beam. 

This time evolution is the subject of this paper. In order to have a clear theoretical model which 
is exactly solvable we restrict to the case of a two-level atom moving in a running laser wave. Despite 
its relative simplicity the two-level model is well suited for the description of certain experiments such as 
Ramsey spectroscopy or atomic interferometry H, || . But this exact solvable model is not only accessible 
to experimental investigation. It may also serve as a reference for more complicated atoms with more than 
two levels which cannot be solved rigorously. Numerical methods used in this context can be tested with 
the two-level model. Another important advantage is the detailed physical discussion which can (and will) 
be made using exact solutions. The effect of every term in the Hamiltonian can easily be identified. 

We describe the dynamical behaviour of the two-level atom under the simultaneous influence of the 
running laser wave and the earth's acceleration in the Schrodinger picture by means of the unitary time 
evolution operator 

\m) = u(t)\m) ■ (i) 
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which is in our case a 2 x 2 matrix. It is methodically essential for the following that we follow an algebraical 
approach to derive U (t) . This approach makes only use of commutation relations and turns out to be 
particularly transparent as far as a continuous compelling physical interpretation is concerned. Because of 
its independence of the initial state |^>o) the generality of the approach is evident. Moreover, it turns out 
that the resulting exact expressions for U(t) are very compact. We mention that the algebraic method has 
already been used in the context of atomic interferometry to calculate the phase shift induced by external 
potentials including gravity for very general interferometer geometries . The gravitational phase shift was 
also calculated by Borde Q with a similar method. 

To get a clear picture of what is going to happen it may be useful to refer to the equivalence principle. 
It states that the influence of a homogeneous gravitational field on the atom moving in a laser wave can 
be simulated by constant acceleration. This means that the following situation is physically equivalent to 
our original set up: A two-level atom is at rest or moving with constant velocity relative to an inertial 
system. The laboratory with the laser attached to it moves with constant acceleration. The consequence is 
that the laser wave reaches the atom with Doppler shifted frequency. Because of the acceleration this shift 
changes in time. It acts as a time dependent detuning. The internal transitions of the atom in the limit 
of vanishing acceleration and accordingly constant Doppler shift are essentially known. For not too large 
detuning they are described as Rabi oscillations. In our case, because of the time dependent detuning, the 
internal behaviour of the atom in time must be worked out anew. This is the central task of this paper. 
Note that there is no further influence of gravity on the internal dynamics of the atom, because the fact that 
reference is made to an accelerated reference frame can have no measurable consequences in this regard. 

Turning to the center-of-mass motion of the atom we know that in the rotating wave approximation the 
transitions between the two internal states are related to the absorption or emission of one photon. The 
corresponding energy and momentum transfer by recoil must show up in the final state \ip(t)}. It is contained 
in a simple way in the transition matrix elements Uvi{t) and U2i(t). 

Finally, there remains the motion of the atom relative to the accelerated reference frame (or from the 
other point of view, the free fall in the gravitational field). We will separate this center-of-mass motion 
from the internal dynamics and the recoil effects in factorizing U(t) with exp(— iif c . m .i/?i) as a factor on the 
left. H c .m. is the center-of-mass part of the Hamiltonian. Clearly this factor drops out if one calculates the 
probability of finding the atom in the excited state. This probability is not related to the actual position and 
momentum of the atom. It is influenced by gravity and acceleration only via the fact that the atom registers 
the laser wave during the time t with a time dependent Doppler shift. As stated above, the additional fact 
that the description refers to an accelerated reference frame does not influence internal atomic processes. 

The paper is organized as follows. In section 2 we will introduce the model Hamiltonian and perform 
a unitary transformation to get a time-independent Hamiltonian. The exact result for the time evolution 
operator U(t) will be derived in section 3. Sections 4 and 5 are devoted to a small influence of gravity 
and the long time limit, respectively. A surprising and not commonly known breakdown of perturbation 
theory will be demonstrated in section 6. In section 7 we will conclude the main results and discuss possible 
applications. 



2 The model 

We suppose the Hamiltonian of the two-level atom to be of the form 

„_(*-!*/» Ea . iV2 ) + ( i .y^-i.s, (2 , 

where E e (E g ) and 7 e (j g ) are the energy and the decay rate of the excited state \tp e ) and the ground state 
\i/) g ), respectively. Although we have inserted the decay factors in a phenomenological way, a more complete 
description of the spontaneous emission would require a density matrix approach. For this reason the validity 
of our model is restricted to times shorter than the life time of the states, but this may be quite long for 
metastable states. For convenience, we will set 7, = in the calculation. This is no loss of generality since 
they can be reintroduced by replacing Ei by Ei — i"fi/2 in all expressions. 

H c . m . = -^-Ma-x (3) 



2 



is the center-of-mass part of the Hamiltonian with x, p, and M being the position and momentum operator 
and the mass of the atom, a denotes the constant gravitational acceleration acting on the atom. The 
interaction with the running laser wave is modeled by the dipole coupling 



d- E — Ml cos[aj L t - k ■ x + y>] ^ 1 j , (4) 

where O := d- Eq/% is Rabi's frequency. Eq is the amplitude of the laser wave, u>l its frequency, (p its phase, 
and k its wave vector, d is the dipole moment of the two-level atom. 

In this form the Hamiltonian depends explicitly on the time t. To get rid of this dependence thus making 
the dynamics simple we make a unitary transformation with the operator 

0{t)= (eM^Lt-k-x + V ),2\ \ (5) 

V exp|— i\ujjj, — k ■ x + f)/2\ I 

In the rotating wave approximation, i.e., after neglecting all terms oscillating with the frequency 2u)l, we 
find for the transformed Hamiltonian H = OHO^ 1 — ifiOO^ 1 the expression 



E g ) 1 \~ c m - ' 4 J V 1 / 2 
where we have introduced the well known recoil shift 



hk 2 

d := 

The operator 



D:=^p-k (8) 

is crucial for the following calculations and their physical implications. For an atom moving with velocity 
v = p/M it can be written as v ■ k. Hence this operator represents the Doppler shift of the laser frequency 
in the rest frame of the atom. We will call it the Doppler operator. Strictly speaking we are in this context 
only allowed to argue with reference to the transformed momentum operator 

p\ = OplO- 1 =pl + ^hka 3 (9) 

where 1 is the unit matrix in two dimensions and o~i are the Pauli matrices. But with Eq. (^) it is easy to 
see that the term containing the Doppler operator can be written as 

Da 3 = Da 3 - SI (10) 

so that the interpretation is the same apart from a trivial shift of the energy eigenvalues. 

In addition to the fact that it implies a time independent Hamiltonian, the unitary transformation (^) has 
two important consequences. The momentum transfer to the center-of-mass motion related to the absorption 
or emission of a photon has been absorbed in 0(t). That the laser causes internal transitions is reflected 
as usual by the O-term. The Hamiltonian (||) shows very clearly that after separation of the momentum 
transfer the term containing the Doppler operator and therefore the operator p is the only one which couples 
the internal degrees of freedom to the enter-of-mass motion in a nontrivial way. 



3 The exact time evolution operator 

The main advantage of the unitary transformation (^|) was to get rid of the explicit time dependence in the 
Hamiltonian (after the rotating wave approximation). As a consequence the new evolution operator U(t) 
can simply be determined by calculating 

Oitpi^O-^O) = U(t) = exp[-itH/h] =: e A+B (11) 
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with 



A := Ql 

B := Pda + Ecri (12) 



and 



Q := ~Y \ H cm. + \ [E e + E g + US/2] 



P := |{A-^} 

R ■= fn. (13) 

The operators A, B, P, Q, and R are introduced to clarify the mathematical structure of the calculation 
below. In P the quantity 

E P . — E n 

T 

denotes the detuning of the laser frequency with respect to the atomic transition. 

To separate the free fall of the atom's center of mass we factorize U(t) according to 



A := uj l - uj eg with uj eg := — ^ — - (14) 



U(t) = e A W{t) . (15) 

To know the complete time development of the two-level atom it now remains to determine the operator 
W(t). It contains the influence of the gravitation on the internal dynamics. It is easy to recognize its 
structure. The factorization ( [l5| ) will introduce in W(t) the commutator [^4, B] and therefore 

[Q,P] = l-[D,H c . m ] = ~k-at 2 . (16) 

This is a c-number. The second equation shows that for k ■ a ^ gravity causes a time dependent Doppler 
shift. This is the central physical effect. We introduce 

D t := D + k-at (17) 
(it may also be written as k ■ (p/M + at)). The relevant time scale is 

= ■ (18) 



k ■ 



For an optical laser with \k\ « 10 7 m 1 and the earth's acceleration (|a| = 9.81 m/s 2 ) r a is about 10 4 
seconds. Introducing 

C := sgn(£ • a) (19) 

we can rewrite Eq. (|l6|) as 

[Q,P]= l -£. (20) 

z ~a 

To work out W(t) of Eq. ( |l5| ) we employ a method which was used by Lutzky |9| in the context of the 
Baker-Campbell-Hausdorff formula. We consider the operator 

G(X) := exp[A(A + B)] =: e XA W(X) . (21) 

and restrict to A = 1 at the end. Differentiation of Eq. ( pl| ) with respect to A leads to the differential equation 

^ = ( e ~ XA B e XA ) W{\) (22) 
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with the initial condition W(X — 0) — 1. Using the identity 

K B e XA = V — X„ , K := B , X n+1 := [if„, A] (23) 
z — ^ n \ 



-XA D „AA 

n! 

71=0 



which holds for any two operators A, B one can simplify Eq. ( p2| ) in our case to 

dW 



i/x (B-X[Q,P]a 3 )W(X). (24) 



It may be written as a matrix equation, 



dWn/dX dW 12 /dX\ _ f P-X[Q,P] R \ ( W n W u 

dW 2 i/dX dW 22 /dX I ~ { R -P + X[Q,P] I I W 21 W 22 



(25) 



This is an operator-valued system of differential equations, but it contains only commuting operators (since 
[Q, P] is a c-number) so that we can treat it as an ordinary differential equation. 

Inserting the equation for dWn/dX (for dW 22 /dX) into the equation for dW 2 \/dX (for dWn/dX) one 
arrives at 

' hW - = {R 2 -[Q,P] + (X[Q,P]-P) 2 }W 11 



dX 2 
t 2 W 22 



dX 2 

After the introduction of the parameter 

R 2 . n 



= {R 2 + [Q,P] + (X[Q,P]-P) 2 }W 22 . (26) 



(27) 



^ h 4 + 0-~}Wn(y) (28) 



2[Q,P] 4k -a 

and the change to the variable y := (A[Q, P] — P)y/2/[Q, P] Eq. ( |2^ ) becomes 



dy 2 { 4 ' " 2 



dy 2 t 4 ' " ' 2 

with the initial conditions Wu(X = 0) = W 22 (X = 0) = 1 and 

dW u 



yl + 9 + l -\W 22 {y) (29) 



dy 



dW 2 



^ = Wra = ^ TQA0 (30) 



(Aq := A — D), We used the opportunity to introduce the operator of the time dependent Doppler shifted 
detuning 

A t :=A-D t = A-D-h-at (31) 

in which the effect of acceleration is already included. 

The solution of Eq. ( p8| ) and Eq. (^9|)is a linear combination of parabolic cylinder functions. All relations 
between these functions used in the following are taken from chapter 19 of Ref. fTofl . The standard solutions 
are given by 



u( ]= r- -f/i / iFi(a/2 + l/4, 1/2, y 2 /2) _ y ^(q/2 + 3/4, 3/2, y 2 /2) 
[a,V> V?re 1 2«/ 2 +V4r( a /2 + 3/4) 2«/ 2 - 1 /4r(a/2 + l/4) 



(32) 



and 

V(a, y) = r(a + V2) {sin(7ra)[/(a, y) + U(a, -y)} (33) 
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where iFi(a, f3, y) is the confluent hypergeometric function. For W\\ we have to set a = — 1/2 + 9, and for 
W 22 we have a = 1/2 + 9. The linear coefficients can be derived from the initial conditions. By using the 
Wronskian relation 

vf-f-V-Jl (34) 
ay ay \ ir 

as well as 

^(a,y) + ^yU(a,y) + (a+^)U(a + l,y) = 

^-(a,y)+ l -yV{a,y)-V{a + l lV ) = (35) 
one deduces (and from now on we set A = 1) 

Wn(p) = J^{v(h + e,]^r„Ao)C/(-i + 9, ^TaAt) + 0U($ + 9, ^L Ta A )V(-± + 6, ^r a A t )} 

W 12 (p) = J^{u(-i + 0, 7 =r Q Ao)F(-i + 9, ^Lr a A t ) - V(~i + 9, ^L Ta A )U(-i + 9, ^r a A t )} 

W 21 (p) = ^9{u(i+9^T a A )V(i + 9^T a A t )~V(i + 9,^T a A )U(i+9, 7 i.T a A t )} 

W 22 (p) = J|{^(-i+fl J ^r„Ao)^(i+fl J ^r A t ) + ^(-i + fl,^T Ao)ni+«»^aAt)} (36) 

The time dependence is contained in A t . 

With Eq. ( |36"1) the total operator U (t) is given by 

E e + E g _ 5 , 1 TT \\ ( W u (p) W 12 (p) 



U{t) = exp < -« I ^ + - + -H^j I m ~ wJ{p) . (37, 



It remains to cancel the initial unitary transformation 0(t) in Eq. ( |ll| ) to obtain the exact expression for 
the time development operator: 

U(t) = o^i^u^oio) 

E e +E g , 5 , 1 „ M ( e -it(u>L~D)/2 e ik-at 2 /4 



(38) 



6XP + 2 + ft^ J / ^ e «( Wi -iS)/ ae ^S.«-/4 J x 

/ Wu(p-hk/2) Wu(p- hk/2) e^-*) \ 

\W 21 {p + hk/2)e~^- s -v) W 22 {p + M/2) )' 



This is the main result of the paper. The argument p± hk/2 in the operators Wij denotes that the operator 
p has to be replaced by this expression wherever it occurs in of Eq. (|36|). 

For practical calculations it is useful to rewrite the factor exp{— itH c . m Jh} in Eq. ( p8| ) with the aid of 
the Baker-Campbell-Hausdorff formula. Doing so one arrives at 

ex 4^l^} exp {f M ^4 ex 4l^4 exp {^ Ma "} ■ (39) 



This form is especially convenient for the interpretation. The first exponential in Eq. (|39|) describes the 
ordinary kinetic energy of the atom. The second term describes the momentum gain M at of the accelerated 
atom since the operator exp(ifc • x) for any c-number vector k amounts in adding Tik to the momentum. 
The third term is responsible for the displacement of the atom by the amount at 2 /2, and the last term is a 
time dependent c-number phase factor which corrects the "error" that we have made in writing all factors 
of H c . m . in separated exponentials. 



G 



Turning to the discussion of the result ( pq ) we note that the momentum transfer to the atom related to 
the absorption or emission of a photon can be read off from JJxi and XJ 2 \ (see the exponential factor on the 
right). The time dependence of W goes solely back to the time dependent Doppler operator D t contained 
in At which modifies the detuning. D t introduces thereby a dependence from the center-of-mass state of 
the atom via the momentum operator. Those who are familiar with the evolution of a two-level atom in a 
running laser wave may miss the recoil shift 5 in At which normally comes with the detuning [| . It enters 
At only after the retransformation of U in Eq. (|38|). The replacement of p by p ± hk/2 in Wij(j) ± hk/2) 
amounts in replacing At in by A t =p d so that the missing recoil shift is reproduced. Rabi's frequency Q 
appears in the parameter 9 in the form f2r a . It grows with the intensity of the laser wave in relation to the 
characteristic time r a . 

An instructive consistency check is to turn off the laser by setting f2 and hence 9 to zero. In this case 
the freely falling two-level atom should be recovered. 9 = implies immediately W\ 2 = W21 = 0. For the 
evaluation of the Wa it is necessary to use the identity [[To| 

U(-l/2,y) = e-y 2 l i , I/(l/2,y) = /V/ 4 . (40) 



Inserting this into Eq. (36) and combining the result with Eq. (38) leads to 



U{t) = exp <{ —H c . m . \ I . Eat/H ) (41) 



% cm j \ e -iE g t/n 

as was to be expected. The first term describes the free fall of the atom, and the second term contains the 
internal oscillations. By using Eq. ( pp| ) it is also possible to derive an expansion of Eq. ( p8[ ) for small Q. But 
since the result contains various combinations of Error functions and is not much better to interpret as the 
complete result we will omit it here. 

We will close this section with a mathematical note. In the definition of ^fiQ the square root of a complex 
factor appears. Since ^fz for complex z is a multivalued function we should decide which branch has to be 
taken. Fortunately, this not a big problem since the variable y in Eqs. ( |32"1) and (|33j) (and hence in 
Eq. (^)) is squared in the argument of the function There is only one linear factor of y in Eq. ( |32| ) 
which can give an additional sign. 



4 The limit of small gravitational influence 

We now consider the limit of very large r a , or more precisely small [Q, P]. The results will thereby be valid 
for all times t with t <C t q so that the atom will be able to perform many Rabi oscillations. The physical 
meaning of the limit is quite clear. Very large r a means very small k ■ a, i.e., the momentum transfer from 
the laser beam to the atom is only slightly altered by the acceleration because the latter is very small or 
nearly perpendicular to the laser beam. The technical details of this limit are somewhat involved and are 
explained in Appendix A. The resulting expansion of Wij (p) up to linear powers of k ■ a is 



Wnip) 

W X2 (p) 
W 21 {p) 

W 22 (p) 



A 

cos(o>t/2) + i^r- sin(d>f/2) + 

LJ 



-k ■ at 2 



A 2 



in 



cos{£jt/2) 



n 



sin(a)t/2) 
u>t/2 



= — sin(a>t/2) k ■ at-^ I H A 



it ; 



in 



n 



— sin(cDi/2) + -k ■ at— 1 - — A 



it 



cos(o></2) — i^- sin(wt/2) + 



-k-aA-^ 



2 r 



cos(Cjt/2) 



sm(ujt/2) 



jt/2 



A . /A sin(wt/2) 



cos(u)i/2) 
cos(ujt/2) 



LOt/ 2 

sin(wt/2)l 

Cot 1 2 J 
$m{Lbtj2)\ 

cot/ 2 J 



.A . sin(d)t/2) 

LO LOt/ 2 



(42) 
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Here we have defined the frequency operator 



^2 + (-A ) 2 (43) 



which incorporates the well known fact that the frequency of the Rabi oscillations is altered when the laser 
frequency is detuned versus the atomic transition frequency. Note that also the effect of the Doppler shift 
is included in Ao. The additional "-" sign before Ao indicates that the negative branch of the square root 
has to be taken if f2 is set equal to zero. 

The analysis of Eq. ji^ ) is relatively simple. First, it is easy to show that for k ■ a = Eq. ([l||) describes 
the ordinary Rabi oscillations of a two level atom in a running laser wave. To see this we solve Eq. (^4|) with 
[Q, P] = which corresponds to k ■ a = 0. The obvious solution is 

W = exp{AS} . (44) 

Since in the operator B = P03 + Ra\ the operators P and R commute one can apply the formula 

exp{— is ■ (?} = 1 cos(v / i^') — i— = sin(V / i^) (45) 



to reproduce just the a-independent part of Eq. (]42|). This result describes in operator form what has been 
obtained in previous calculations |^] where the time evolution was derived (for certain wave packets) in 
momentum space. It is interesting to see the effect of the Doppler shift on the atomic evolution. It not only 
alters the frequency uj of the Rabi floppings. A large Doppler shift also damps the transition probability 
because of the l/u) dependence of the matrix elements W^ip) and W2i(p), and it is (together with the 
detuning A) responsible for the imaginary part of W\x and W^i- 

Since in our model gravity alone cannot cause transitions, the corrections to the transition matrix elements 
W12 and W21 have to vanish if the influence of the laser disappears, and this is exactly what happens in 



Eq. (42). We have seen this already above. It may be surprising why we have in this limit corrections at 
all, because for the almost freely falling atom the influence of the homogeneous gravitational field is already 
included in the factor of exp{— itH c , m ./h} in Eq. (|38|). The remaining corrections are indeed only necessary 
to cancel the phase factors of exp{±ik ■ at 2 / 4} appearing in Eq. (p^). 

Like the unperturbed part the corrections for small k ■ a are oscillating with the frequency uj/2. Further- 
more, it is not difficult to see that all corrections are bounded functions of the operator uj. This enables us to 
study their time evolution. For t — they vanish. For small t the corrections to the diagonal matrix elements 
of W grow like t 2 , and the non-diagonal corrections grow with t 3 . For moderate t this cubic dependence 
becomes a linear envelope of an oscillating function. The corrections are all suppressed if the Rabi frequency 
becomes large. In this case the force caused by the laser beam is much larger than the gravitational force 
so that the former dominates the evolution. Note that the corrections in W11 and W22 do not vanish for 
large Doppler shifted detunings. This is a reasonable result since a large detuning means that the laser is out 
of resonance and does not affect the atomic evolution anymore. At the same time the influence of gravity is 
not altered so that there must remain some effect of gravity even when the detuning is large. 

Another effect of gravity can be anticipated by examining W 12 . For vanishing corrections this is a purely 
imaginary operator. Switching gravity on we see that it develops a real part which can become large if t 
grows. This is an indication that the acceleration causes additional time dependent phase factors. We will 
look at this more closely in the next section. 



5 The evolution for long times 

The results of the previous section have been obtained in discussing the weak gravity limit of the exact result. 
They could have also been worked out within an approximation scheme which starts with the gravitation- 
free case and perturbs it by a small k ■ a. In this section we focus on a result which is of non-perturbative 
character: the behaviour of the atom at late times t. The formal condition which we will use is 

|A t | = \A -k-at\ » — . (46) 
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In this case we can apply the asymptotic form of the functions U and V in Eq. ( pq ) which depend on A*. 
Since A t is an operator it would be more precise to require |(A t )| 3> l/r a and to consider only wavepackets 
for which the standard deviation of A t remains sufficiently small. Because any time dependence of Wij is 
enclosed in A t it is easy to see that the condition ( [46| ) is equivalent to f > t„ (if Ao is not too large). 
Physically this means that we consider the case when the Dopplcr shifted detuning becomes large if it was 
initially (this is Ao) not very large. It would be a different physical situation if the detuning is initially very 
large and becomes small due to the Doppler effect of the accelerated atom. This will not be considered here. 

The main ingredients of the limiting process in question are explained in Appendix B. The matrix W is 
found to be 

. iCU(i + 6, -^T a A ) 

Wn(p) « \/ ^^(- C r a & t) -WJ* e ^n° \ V (i+e, ^Lt„ A ) - f( -^ 

n„(|7i * \/^e^ A */ 4 (-Cr a AO^^/S-^/ 16 C/(i + ^^T a Ao) 
4fc • a v 

W 22 (p) « e^ A ?/ 4 (-Cr a A t ) 4Co2 ^ 2/4 e -"° 2T » /16 t/(^ + 0, ^r.Ao) . (47) 

The most striking feature of these operators is that (to lowest order) only their phase varies with A t and 
hence with t. This is not exactly true because terms containing A t are operator valued and can alter also the 
shape of a wavepacket. But for sufficiently narrow wavepackets in momentum space the factors containing 
A t in Eq. ( f47j ) simply produce an additional phase shift with logarithmic time dependence (remember that 
y lc — exp(iclny) ). From Eq. (|38| ) we then can conclude that with growing time the matrix elements of the 
total evolution operator U (t) vary also only in their phases. We will discuss this time dependence first and 
turn to the amplitudes thereafter. 

We interpret the result as follows: Atoms exposed to gravity and a laser beam may have been in resonance 
with the laser at some earlier stage of their evolution. In this stage they perform a number of Rabi oscillations. 
But since the momentum dissipation due to spontaneous emission is neglected in our model, the atoms loose 
during each Rabi cycle the same amount of momentum as they gain, the net effect being zero. During the 
Rabi oscillations the atoms are accelerated by the earth's gravity so that their velocity increases. Because 
of the Doppler effect the atoms are then driven out of resonance with the laser beam so that the Rabi 
oscillations are vanishing with increasing time. The transition between excited and ground state is frozen. 
This is reflected by the fact that the absolute value of the is almost constant for long times. 

Looking more closely to the particular matrix element U\\ of Eq. (p^) we can write the time dependent 
phase factor as 

CXP {"*' (^T* + \ B - + Q } CXP {^P 1 + a 3 } e-^^^) 2 / 4 [- C r a( A t+5) ]-^^/4 . 

(48) 

Expanding the factor of A t and keeping only time dependent terms leads us to 

cxp {-^f 1 } exp j-ftiff c . w .} [-Cr„(A t + S)}^^ 4 . (49) 

Here we see that each phase factor (but the last) is linear in t. The time dependence of U12 agrees in this 
limit with the one of U\\. This is reasonable because both matrix elements correspond to atoms which are 
excited at the time t. The result for U21 and t/22 can be read off from Eq. ( [49| ) by replacing E e by E g (these 
matrix elements describe atoms which are in the ground state at time t) and [ — Qt a (At + S)]~ i(n r » /4 by 
[-Cr Q (A t -S)} 1 ^ 2 ^/*. Each factor in Eq. © allows a proper physical interpretation. The first exponential 
describes just the internal energy of the excited state. The second is the time evolution of a free point 
particle in a homogeneous gravitational field. Hence we see that the atom is essentially freely falling if the 
third term is neglected. 
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Nevertheless, this slowly (logarithmically) varying phase is of physical interest since it contains both 
and a and is therefore the only remaining time dependent effect arising from both gravity and laser light. 
Because it contains Q its origin must be the fi-term in the Hamiltonian (Q) which induces transitions between 
ground and excited state. But since it is only a phase factor it does not describe such transitions (which 
would be characterized by a population transfer between \ip e ) and , i.e., in a time dependent change of the 
absolute values of W%j). Seen in this way the result seems to contradict our intuition, but this contradiction 
can be resolved by a comparison of the present situation with Raman transitions (see, e.g., Ref. pi). Raman 
transitions are possible in a system with two lower and one upper state ("A system"). In this A system one 
can induce direct transitions between the two lower states without populating the upper state by applying a 
laser field with a large detuning versus the transition frequency between a lower and the upper state. This 
is similar to our case: we also have a large detuning and no population transfer to the upper state, the only 
difference is that we have as a form of degeneration only one lower state. Hence, we interpret the third 
term in Eq. (^9|) as a kind of Raman transition in a two-level system for which the large detuning prevents 
the population of the upper state during a Rabi cycle. The peculiar logarithmic time dependence of the 
phase is a consequence of the time dependent detuning. We should mention that because we have neglected 
spontaneous emission the argument does not only apply to Rabi cycles of the form \ip g ) — > \ip e ) — > \ip g ) but 
also to cycles of the form \ip e ) — > \ip g ) — > \ip e ). 

We turn now to the discussion of the time independent part of . The phase of this part is a rapidly 
oscillating function of Ao. To get a feeling for the absolute values we consider the case of a narrow wavepacket 
of atoms with an initial velocity so that Ao vanishes, i.e., the Doppler shifted detuning is initially zero. With 
the aid of the formulae in appendix B (especially Eq. (|6§|)) it is not difficult to obtain 

|Wu(Ao=0)|= |W 22 (A = 0)| = ./!(! + e-^ 2 /4) 

V (50) 

|VMA o =0)|= |W ai (Ao=0)|= ^(l-e-^/4) 

This is a surprisingly simple expression with several interesting features. First, it only depends on the 
absolute value of k ■ a since the variable £ is absent. This is reasonable since the appearance of C^ 2r o m 
the exponent would allow imaginary absolute values of W\% and W 2 \ ■ Second, for vanishing laser intensity 
(p, = 0) we find \Wiz\ = |H^i| = and \Wu\ = |W / 22| = 1 in accordance with the fact that without laser 
beam any internal transition is impossible. For high laser intensity (fl 1) we get a complete mixing: 
| Wij | = 1/ \/2. This can be understood as a consequence of the large frequency of the Rabi cycles. For 
Q — > oo there is an infinite number of Rabi cycles per unit time. This is not a well defined result. But 
one can replace these quick oscillations by there average, and this procedure results of course in an equal 
probability for ground and excited state as indicated by |Wij| 2 = 1/2. In general the transition amplitude 
depends on the variable fi 2 r 2 which can be interpreted as representing the influence of the laser (because of 
Q) acting for the time r a after which the atoms are essentially out of resonance. 



6 The breakdown of the Magnus expansion 

We now have finished the physical discussion of the atomic dynamics in a homogeneous gravitational field 
and a running laser wave. This last section of the paper is of purely theoretical interest. We will show 
here that the application of the Magnus perturbation expansion jL2| would lead to unphysical results for the 
evolution operator. This fact was already examined in the literature (see Ref. [|l3| and references therein), 
mostly for the harmonic oscillators and simple two-level systems. A further analysis for the falling two-level 
atom in a running laser wave has the advantage that it is possible to use the Schr" odinger picture and that 
we can compare the result of the Magnus expansion with the exact solution obtained above. 

To perform the Magnus expansion we go back to Eq. ( p4| ) and ask whether it is possible to treat [Q, P] 
as a small term and to apply a perturbation expansion instead of solving Eq. (p4j) exactly. The Magnus 
expansion consists of an expansion of the exponential of an operator. Setting W — exp(.F(A)) we may try 
to calculate F to first order in [Q,P]. This can be done by using the equation (see, e.g., Ref. Q) 

A °° ( -n(J+i) 
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with K\ :— dF/dX and Ki + \ := [Ki,F] which is exactly valid for any operator F provided the exponential 
makes sense. Eq. (|5l]) has the same structure as Eq. (Q). The strategy to solve Eq. ( p4| ) perturbatively is 
therefore to find an operator F such that the sum in the r.h.s. of Eq. ( |5l| ) reproduces, to first order in [Q, P], 
just the prefactor of W in the r.h.s. of Eq. (p4|). In view of this prefactor we make the following guess for F: 

F = XB- -X 2 [Q, P](7 3 + [Q, P] J2 , ' b nRn (52) 

2 ^ (n + 1)! 

where b n are numbers, R± := [<73,P], and R n +i '■= [R n ,B]. Inserting this in Eq. ( |5l| ) and comparing the 
result (to first order in [Q,P]) with Eq. (E4J) leads to the condition 



n— 2 / \ 

(n + 2)6„ - - + ( ™+ 2 ) bn-1-1 = , n > 2. (53) 

It is easy to check that the numbers b n are related to the Bernoulli numbers B n by b n = B n+ \ by using the 
relation 

£(?U = 0. (54) 

A closed expression for F can be found by noting that 

R 2n = (itui) 2n - 2 R 2 , R 2n+1 = {ttiu) 2n R 1 (55) 
and i?2«+i = hold. Putting everything together one finds 

F(X) = XB - y[Q, P]a 3 + [Q, P]Ri^ {J£^J - 1 + ^} + 0{[Q, P] 2 ) . (56) 

The approximate solution of our problem is obtained if A is set to 1. It is obvious that this solution becomes 
singular whenever the condition 

tot = 2irN , N e N (57) 

is fulfilled (this is possible if we neglect the decay rates j e and j g ; if they are non-zero the solution P(A) 
contains an unphysical resonance). 

The approximation should be good for \[Q, P]\ <§; 1, i.e., for t r a . But u> is for small detunings 
essentially identical to Rabi's frequency fi. Hence the first singularity occurs at t = 2ir/fl which is much less 
than r a for not too small Q. This illustrates that a perturbation approach based on the Magnus expansion 
fails to describe our system properly. This fact was explained as a consequence of a finite convergence radius 
of the Magnus expansion [ p3| . It may be that the reason lies in the splitting of the right hand side of 
Eq. ( pl| ) into one factor e F which contains any power of [Q, P], and in the sum over I which in our example 
is calculated only to first order in [Q, P]. 



7 Conclusions and Outlook 

In this paper we have exactly determined and discussed the time evolution of a two-level atom falling in a 
homogeneous gravitational field under the influence of a running laser wave. The time evolution operator has 
been worked out in an algebraical way. For neglected spontaneous emission the Doppler effect was shown to 
be the origin of the characteristic new physical features. An atom which is initially in resonance with the laser 
beam will first perform Rabi oscillations. The homogeneous gravitational field accelerates simultaneously 
the center-of-mass motion. During this acceleration the Doppler shift causes the laser frequency to drive out 
of resonance with the atom. Thus the Rabi oscillations are fading away until only a slow variation of the 
phase remains, no population transfer appears for large times. This situation is similar to Raman transitions 
in A systems. 

The treatment of the model described above may be regarded as a first step towards the inclusion of the 
influence of gravity in situations of greater practical importance. That there is (almost) no net momentum 
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transfer of the laser to the atom relies on the neglection of the spontaneous emission. As long as the laser 
is in resonance with the atom each Rabi cycle will be completed, no total momentum transfer occurs. But 
if spontaneous emission is included Rabi oscillations can be incomplete. The resulting momentum transfer 
from the laser to the atom may be exploited to construct a gravitational atom storage if it is adjusted so that 
it can cancel the acceleration by the homogeneous gravitational field. The interruption of the Rabi cycles 
can also be managed in a coherent way, e.g., by using a three level "A" scheme with two hyperfine ground 
states. In this case a magnetic field (see, e.g., Ref. |15|) or a microwave may be used to carry the atoms back 
in their initial state. It should be mentioned that there is already a gravitational atom trap |l6| which was 
theoretically studied by Wallis, Dalibard, and Cohen- Tannoudji |l7]]. The difference to our proposal is that 
we use the simultaneous action of gravity and laser forces whereas in the existing device these forces act in 
different time periods. 
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A 

The aim of this appendix is to study the behaviour of U(9 — 1/2, yo + e) for small 

e := v/^QTpI ~ - (58) 

and to derive the corresponding limiting case of Wij. For brevity we have introduced the parameter yo ■= 
—iT a Ao/i/i(. We start with the observation that any derivative of U(8 — 1/2, y) with respect to y can be 
written in the form 

(9 ~ 1/2, y) = X n (0 - 1/2, y)U{6 - 1/2, y) + Y n (0- 1/2, y)U'(9 - 1/2, y) (59) 
where the prime denotes the derivative with respect to y and the functions X n and Y n fulfill 

Xn+i = X' n +(^+e-l/2\Y n 

Y n+1 = Y^ + X n (60) 

with Xo = 1 and Yq — 0. This relation, which follows directly from the differential equations (p8| ) and 
(p9|), will turn out to be useful for the Taylor expansion of 

oo i 

U(9 - 1/2, y + e) = £ (<? - 1/2, Vo) (61) 

around yo- F° r the expansion of X n and Y n in e it is of importance to note that 

— 2P R 2 

y = — , e = \ (62) 

holds. Bearing this in mind one can use Eq. ( |60| ) to proof by induction that 

_-2 



X 



2/i 



2 "{(itd>/2) 2 " T ^e 2 (itoj/2) 2n - 2 + 0(e 4 )} 



Y 2n = -Pe- 2 "+Vn-l) j(^/2r"~ 4 T- ^-{itu/2Y n - b + <9(e 4 ) 

A 2n+1 = -Pe- 2n+1 n 2 ((z^/2) 2 "- 2 T ( " ^/2) 2 "' 4 + 0(e 4 ) 



Y 2ll + I = : \(itu/2) 2n *^(itu/2) 2n - 2 + 0(e 4 )} . «(«! 
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is valid. Here we used the fact that 



,2 



l -tu> = ssJe + v -± = ^(-py+w (64) 

is independent of e. The operator Co is defined in Eq. (fi"3|), and the upper (lower) sign in Eq. (63) holds for 
U(6 — 1/2, yo) and U(6 + 1/2, yo), respectively. Note that the same expansion holds for V(6 — 1/2, yo) and 
V(9 + 1/2, y ) simply because they are also solutions of Eqs. ( |2g| ) and ( p9h, r espectively, and because these 
differential equations are the only relation which were used to derive Eq. (|60|). 

We are now ready to calculate the expansion of Wi2(p) in terms of e. We have with Eq. ( |§6| ) 

W 12 (p) = j {U (6 - 1/2, y )V (6 - 1/2, y + e)-V (9- 1/2, y )U (6- 1/2, y + e)} 



£ n-l 



2 ^ n\ 

n=0 



Y n (9 - 1/2, y ){U{6 - 1/2, y )V'(6 - 1/2, y ) - U'{9 - 1/2, y )V(8 - 1/2, y )} 



r y„(^-l/2 ) j/ ) (65) 



n 

71=0 



where we have used Eq. ( p9[ ) for the first step and the Wronskian relatio n (p4[ ) for the second step. It is now 
straightforward to derive a closed expression for W12 by exploiting Eq. (|63[) and the Taylor series of cos and 
sin. The result is given in Eq. (|42]). 

The derivation of the expansion for W21 is similar to the previous calculations and will not be reproduced 
here. For W\\ and W22 one further step has to be included since the functions U and V occurring in the 
corresponding expressions of Eq. (^6|) contain both the parameter 6—1/2 and 6 + 1/2. It is therefore 
necessary to make use of Eq. (|3^) and to deduce expressions for the Wu where either only 6 — 1/2 or only 
6 + 1/2 occurs. After this has been done the calculations are the same as for W12. Inserting all definitions 
one arrives at Eq. (^2|). 



B 

In this appendix we will sketch the derivation of the long time behaviour for the operators Wij (p) in Eq. ( |36| ) . 
This can be handled by using (Eq. 13.5.1 of Ref. @) 

^(a, 0, z) « r(/3) ± mz - a + M e ^«-^ (66) 
I (p — a) 1 (a) 

where z is a complex number with \z\ 3> 1, and the upper sign holds for —tt/2 < arg z < 3ir/2 whereas the 
lower sign holds for — 37r/2 < arg z < —it/2. We apply this formula to the function U{— 1/2 + 6, — ir a /\ t / VK) 
and in the same manner to any other function in Eq. ( |36| ) which depends on At by inserting it into Eqs. (^) 
and (^3|). Note that the limit depends on £ since f = ±1 determines arg z. We assume (— r a At) to be a 
large positive number and choose arg — C^/2. Hence the "+" sign in Eq. ( |66| ) holds for C = 1 and the 
»_» s jg n for C = — 1. 

The rest of the calculation is straightforward but long. It is useful to apply zT(z) = T(z+1), T(l/2) = y/ir, 
and 

T(ix)T(— ix) — \Y(ix)\ 2 — 7r/[x sinh(7ra;)] 

T(ix + l/2)T(-ix + l/2) = \T(ix + l/2)\ 2 = tt/ coslnW) 

T(2z) = T(z)T{z + l/2)2 2z - 1 / 2 /v / 2^ (67) 

(for real x, see chapter 6 of Ref. |U^) to handle the various factors of the T function arising in the derivation. 
For the case of vanishing detuning the following equations are of use: 

^"'^ = 2 Q / 2 + 1 /4r(3/4 + a/2) 
v , n ^ 2"/ 2 +V 4 sin[ 7 r( Q + l/2)/2] 

y(a ' 0) = r(3/4-a/2) ' (68) 
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